arXiv:astro-ph/0503609vl 29 Mar 2005 


Inhomogeneous models of interacting dark matter and dark energy. 
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We derive and analyze a class of spherically symmetric cosmological models whose source is 
an interactive mixture of inhomogeneous cold dark matter (DM) and a generic homogeneous dark 
energy (DE) fluid. If the DE fluid corresponds to a quintessense scalar field, the interaction term can 
be associated with a well motivated non-minimal coupling to the DM component. By constructing a 
suitable volume average of the DM component we obtain a Friedman evolution equation relating this 
average density with an average Hubble scalar, with the DE component playing the role of a repulsive 
and time-dependent A term. Once we select an “equation of state” linking the energy density 
(fi.) and pressure {p) of the DE fluid, as well as a free function governing the radial dependence, 
the models become fully determinate and can be applied to known specific DE sources, such as 
quintessense scalar fields or tachyonic fluids. Considering the simple equation of state p = (7 — 1) /r 
with 0 < 7 < 2/3, we show that the free parameters and boundary conditions can be selected 
for an adequate description of a local DM overdensity evolving in a suitable cosmic background 
that accurately fits current observational data. While a DE dominated scenario emerges in the 
asymptotic future, with total fl and q tending respectively to 1 and -1/2 for all cosmic observers, the 
effects of inhomogeneity and anisotropy yield different local behavior and evolution rates for these 
parameters in the local overdense region. We suggest that the models presented can be directly 
applied to explore the effects of various DE formalisms on local DM cosmological inhomogeneities. 

PACS numbers: 


I. INTRODUCTION 


Observational data on Type la supernovae strongly 
suggests that the universe is expanding at an acceler¬ 
ated rateP, |2|- This effect has lead to the widespread 
assumption that the inventary of cosmic matter-energy 
could contain, besides baryons, photons, neutrinos and 
cold dark matter (DM) , an extra contribution gener- 
ically known as “dark energy” (DE), whose kinematic 
effect could be equivalent to that of a fluid with neg¬ 
ative pressure. While the large scale dynamics of the 
main cosmic sources (DE and DM) is more or less under¬ 
stood, their fundamental physical nature is still a mat¬ 
ter for debate, thus various physical explanations have 
been suggested. Cold DM is usually conceived as a col¬ 
lisionless gas of supersymmetric particles (neutralinos), 
while DE can be modelled as a “cosmological constant”, 
quintessense scalar fields, tachyonic fluids, generalized 
forces, etc0,0. The standard approach is mostly to con¬ 
sider a Friedman-Lemaitre-Robertson-Walker (FLRW) 
metric, with linear pertubations, making also the sim¬ 
plest assumption that DE only interacts gravitationally 
with DM. However, there are still some unresolved is¬ 
sues, such as the so-called “coincidence problem”, con¬ 
cerning the odd apparent fact that the critical densities 
of DM and DE approximately coincide in our cosmic 
era nil. Aiming at a solution to this problem and bear¬ 
ing in mind our ignorance on the fundamental physics of 
DM and DE, various models have been proposed recently 
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which include asorted forms of interaction between these 
sources [tLH. 1^ ITI| . 

It is customarily assumed that DE dominates large 
scale cosmic dynamics, so that DM inhomogeneities in 
galactic clusters and superclusters can be considered a 
local effect or can be treated by means of linear per¬ 
turbations in a FLRW background. Thus, a reasonable 
generalization of existing models could be to assume in¬ 
homogeneous DM interacting with homogenous DE, so 
that large scale dynamics is governed by the latter. We 
propose in this paper a class of analytic models which 
provide a reasonable description of inhomogeneous DM 
interacting with a generic homogeneous DE source. The 
models are based on the spherically symmetric subcase 
of the Szafron-Szekeres exact solucions of Einstein’s field 
equations for a perfect fluid source^. However, the un¬ 
derlying geometry of the models we present can be easily 
generalized to include non-spherical symmetries or even 
the case without any isometry, since, in general, Szafron- 
Szekeres solutions do not admit Killing vectors. 

The prefect fluid in Szafron-Szekeres solutions in their 
original conception is characterized by a non-rotating 
geodesic 4-velocity field, so that in the comoving frame 
matter-energy density is inhomogeneous, while pressure 
depends only on cosmic time. As we show in section H, 
it is straightforward to re-interpret this fluid as a mix¬ 
ture of an inhomogeneous dust component plus a homo¬ 
geneous fluid. Such mixture have been considered previ¬ 
ously Ellll but in the context of mixtures of baryons 
and radiation. We consider in this paper only the type 
of models examined in El, by assuming the homoge¬ 
neous fluid to describe a generic DE source, while the 
dust component corresponds to inhomogeneous DM, all 
;osrrQl(!aHitlicJ3das=if, reasonable assumption since the dynam- 
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ical effects of quintessence mostly become dominant in 
very large scales, larger than the “homogeneity scale” 
(100-300 Mpc), while DM (galactic clusters and super¬ 
clusters) is very inhomogeneous at scales of this mag¬ 
nitude and smaller. By conveniently rescaling the free 
parameters and calculating relevant physical and geo¬ 
metric quantities, we show in sections III, IV and V 
that the dynamical equation characterizing the models 
is analogous to a Friedman equation in which the av¬ 
erage DM density evolves in the presence of a generic 
(still undetermined) time-dependent cosmological con¬ 
stant or “Lambda field”. By looking at the regularity 
conditions, singularities and asymptotic behavior along 
timelike and spacelike directions, we show in section VI 
that the DE and DM mixture components behave as 
needed for a reasonable cosmological model complying 
with observations: for all fundamental observers the DE 
source dominates over DM in the asymptotic timelike fu¬ 
ture, while boundary conditions determine how the local 
ratio of DE and DE changes along the rest frames of the 
fundamental observers. However, these frames (hyper¬ 
surfaces of constant cosmic time) have zero curvature, 
hence the inhomogeneities in the models are more suit¬ 
able for a study of large scale inhomogeneities than local 
structure formation, since even overdense regions homo- 
geneize and isotropize towards an asymptotic DE domi¬ 
nated scenario. The main observational parameters are 
appropriately defined and calculated for an inhomoge¬ 
neous and anisotropic spacetime in section VII. 

In order to determine the time evolution of the sources, 
we need to assume a physical model, or “equation of 
state” for the generic DE source (the homogeneous fluid). 
Thus, we assume in section VIII a simple “gamma law” 
equation of state of the form p = (7 — 1 ) /i, where p, p 
are the pressure and matter-energy density of the DE 
source. Such an equation of state leads to a DE homoge¬ 
neous fluid evolving like a FLRW fluid with flat spacelike 
sections with a scaling law of the form p oc which 
is compatible with a scalar field with an exponential po¬ 
tential M Although this is a very simple type of DE 
source, it yields analytic forms for the DM density, ob¬ 
servational parameters and allother relevant physical and 
geometric quantities. 

The assumption of a gamma law equation of state fully 
determines the time dependence of all relevant quantities, 
but the free parameters governing spacial dependence are 
specified in section IX. Suitable boundary conditions can 
always be selected allowing for a description of a local DM 
overdense region in a DE dominated cosmic background 
that accurately complies with observational constraints 
on observational parameters: H for DM and for DE and 
the deceleration parameter q. We provide in this sec¬ 
tion a full graphical illustration of the interplay between 
“local” and “cosmic background” effects on these obser¬ 
vational parameters: for example, anisotropy emerges in 
the local dependence of these quantities on the “off-center 
observation angle” ip, while inhomogeneity leads to local 
conditions in the overdense region (DM dominates over 


DE and q is positive) that are different from those of the 
cosmic background: DE dominates and q < 0, as required 
by an “accelerated” universe whose large scale dynamics 
is dominated by a repulsive force associated with DE. 

The issue of the interaction between DE and DM is 
dealt with in section X. We show that the individual 
momentum-energy tensors for DM and DE are not inde¬ 
pendently conserved, thus the models are incompatible 
with these components interacting only gravitationally. 
However, if we assume the DE fluid to be a scalar field 
quintessense type of source, then the models can accomo¬ 
date various prescriptions for a DE-DM interaction, like 
those proposed in the literature niimnn Einally, 
in section XI we present a discussion and summary of our 
results. 


II. FIELD EQUATIONS FOR A CLASS OF 
INHOMOGENEOUS COSMOLOGIES. 


We start with a spherically symmetric inhomogeneous 
spacetime, the spherical subcase of the Szafron-Szekeres 
solutions [IJ, which can be described by the Lemaitre- 
Tolman-Bondi (LTB) metric element 

ds^ = —(? dt^ + ^ dr"^ + dO^, (1) 

\ — K 

with K = K{r), Y = Y{t,r), dD^ = d0^ + s\v? 9 
and a prime denotes derivative with respect to r. As 
matter source we consider a perfect fluid 

= (e+p)n“u''+pg“^ (2) 


where e and p are the matter-energy density and total 
(effective) pressure. In a comoving representation with 
4-velocity u°‘ = 5“ the field equations for m and are 


Gl = Gl = Kp 
—G\ = Ke 


y2 y2 y ’ 


K' YY' Y Y' 
~2YY' ~ YY' ~ y ~ W’ 
K ^ 2Y'Y' K' 

y 2 y 2 YY' YY' ’ 


(3b) 

(3c) 


where k = SttG/c^ and a dot denotes derivative with 
respect to ct. Since the left hand sides of (TTall and (I3bll 
are identical, we obtain the condition Gg — G^. = 0 , which 
yields 


K Y^ 2Y 

Y^ + Y^ + — = ( 4 ) 

an expected result since the comoving 4-velocity for m 
is a timelike geodesic vector and so is incompatible with 
spacelike gradients of p. If we keep p{t) arbitrary, then 
we cannot obtain a simple integral of 0 , but we can 
transform this condition into 

X + -KpX + = 0, (5) 
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where X = while matter-energy density in 

becomes 

4 XX' K 3 K' 

“ 3 XX'X4/3 2 X1/3X'' ^ 

Hence, once we choose p and K, we can find X by in¬ 
tegrating and then e becomes determined by means 
of The solutions of Einstein’s field equations asso¬ 
ciated with O, 0 and 0 are the spherically symmet¬ 
ric subcase of “class 1” Szekeres-Szafron solutions M- 
However, for K ^ 0 there are no analytic solutions for 
the nonlinear equation 0 , while finding a meaningful 
equation of state relating p(t) and e{t,r) is difficult. 


III. A MIXTURE OF DARK MATTER AND 
DARK ENERGY. 


An alternative approach is to consider only the case 
K = 0 (hypersurfaces t = constant have zero curvature) 
and to assume that total matter-energy density decom¬ 
poses as 

e{t,r) = p{t,r)c^ + (7) 


so that we can re-interpret the momentum-energy tensor 
0 as the mixture 

Tab ^ with: 

= pc^u^u\ = {p+p)u^u^+p 

( 8 ) 


So, it is now more natural to choose p and p, as matter- 
energy density and total pressure of a homogeneous fluid 
describing a specific physical system. Then, integrating 
the linear equation 

X+‘^KpX = Q, (9) 


we can determine the metric and the rest-mass density 
of the inhomogeneous dust component 


K, pc 


2 _ ^XX' - KpXX' 

~ Yx' 


( 10 ) 


This approach has been used in the past for modeling 
mixtures of baryons and radiation [l3l Il4| . but it can 
be equally useful to study the interaction between cold 
DM (the inhomogeneous dust: and a generic (yet 

unspecified) type of DE (the homogeneous fluid 
with negative pressure). Because of their construction, 
the two energy-momentum tensors are not separately 
conserved: = —Tyzpb ^ Oj hence we must have a 

non-minimal coupling between DM and DE (see section 
X). 


IV. DETERMINATION OF FREE 
PARAMETERS 


Since 0 is a second order lineal partial differential 
equation, its solutions have the form 

X = V(t) a{r)+W(t) P{r), 


where a(r), f3{r) are the two “integration constants” that 
follow from the integration of 0 and the functions V, W 
satisfy: V jV = WjW = (3/4) Kp. Since we can arbitrar¬ 
ily relabel the radial coordinate r, no loss of generality is 
involved if we relabel these functions as 


V = 

a = 


W = 

f3 = r^/^f 


( 12 a) 

( 12 b) 


where R{t), T{t) and /(r) are arbitrary functions. In¬ 
serting HU and m in 0 and demanding consistency 
with 0 - 0 , we obtain the equations that (given a choice 
of p = p(p)) determine R and T 


np = 

2R 

~ “r ’ 

(13) 

up = 

3R2 

R 2 ’ 

(14) 

f = 

° 

(15) 


where cq is a dimensionless constant, Rq = R{to) and 
Hq can be identified with the Hubble timescale at to- 
The form of equations OS and is identical to the 
field equations of a FLRW spacetime with flat space sec¬ 
tions, this suggest that we identify p and p with variables 
somehow associated with a FLRW background. Once we 
choose an “equation of state” p = p{p) corresponding to 
a specific DE model for the homogeneous fluid (for ex¬ 
ample, a scalar field), we can find R by integrating 
and m, and then T by integrating 113. 

An alternative approach follow by defining what would 
be a “Hubble factor” for an homogeneous fluid: 

H=^, (16) 


allowing us to combine equations OS and d into 

H = -|(p-hp), np = 3i^^ (17) 

so that, once the “equation of state” p = p{p) is pre¬ 
scribed, we get p = p{H). The functions H and R follow 
by integration of d and which is equivalent to 

integrating d and d- 

Once we select / (after solving 1^3^03 oi' iHTIl - llTCll l. 
the metric 0 and dust density d, as well as any other 
geometric or physical variables, become fully determined. 
Under the parametrization given by (IlH and 11211 . the 
metric 0 and DM mass-energy density d take the 
form 


( 11 ) 
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ds 2 


-c^ dt^ + 


{1 + FTf 
(l + /T)2/3 


dr^ +r^ {I + fT f/^ dfl^ 


(18) 


A{3TH + T)fF + 6H{f + F) - 
3{1 + fT){l + FT) 

^ 2coHoRl i[{2TfF + f + F)R^] 

3i?6 1^1 + fT){l + FT) ^ ’ 

where we used 113 and 


satisfying 

2M = -KpY^Y, (27a) 

2M' = neY^Y'. (27b) 

For spherical symmetry and in the appropriate limit, M 
coincides with the ADM and Hawking masses. In par¬ 
ticular, for the metric o in the form of 03 and matter 
energy density o, M can be found by integrating (1^ 


F = f+^rf'. (20) 

Other important quantities are the expansion kinematic 
scalar Q = u°' -a and the traceles symmetric shear tensor 

^ab — ^[a;fc] (^/^) ^ab 


2M = K 


{pc^ + fi)Y^Y'dr = 2M+'^Y^, (28) 

O 


with 


2M 


K C 


pY^Y' dr 


0 = 3H + Zf, 
a°'b = diag[0,-2E, S, E], 

where 

^ ^ / + f+ 2 /fr 

(l + fT) (l + FT)’ 
^ ^ {F-f)t 

{l + fT){l + FT)' 


( 21 ) 

( 22 ) 


(23) 

(24) 


The metric 03 looks like a FLRW line element modi¬ 
fied by the terms containing T, F and /. In fact, all r- 
dependent variables derived above reduce to their FLRW 
forms: p = a = 0, Q/3 = H, if these “pertubations” 
vanish, i.e if either T = 0 or f = f — 0. This homoge¬ 
neous subcase is a FLRW spacetime whose source is the 
DE perfect fluid with matter-energy density and pressure 
given by p and p. In a sense, if /T <C 1 and FT «C 1 
the models would correspond formally to specific exact 
perturbations of FLRW cosmologies. 


= ^R^r^ff[3H{l + fT) + ft], (29) 

where we have assumed M (t, 0) = 0 and used Y and 
Y' given by m, as well as ira and (EHll expressed 
as F = 2{r^/'^fY/[3y/r). Equation it^ yields then the 
Friedman-type evolution equation 

r2 = —= —+ (30) 

in which M. and up respectively play the roles of “efective 
mass” of the inhomogeneous dust (DM) and a homoge¬ 
neous time dependent A term (DE). The mass function 
M can be related to a volume average of the total energy 
density by 


K {pc^ + p) 


K f (pc^ + p)Y^Y' dr 
/ Y^Y' dr 

6M 


6M 

(31) 


Hence, we can identify the averaged DM and DE densities 


V. A TIME DEPENDENT A-FIELD 


It is possible to interpret the homogeneous DE fluid 
as a time dependent “cosmological constant”. Eor this 
purpose we recall that m and m are particular cases 
of the general spherically symmetric spacetime 

ds^ = r) c^dt^ + B^[t, r) dr^ + Y^[t, r) dilY. (25) 

If the source of (ESI is also a perfect fluid like @ , we can 
define a “mass function” M 




^/2 


(26) 


KC^{p) 


k{p) 


6M 41311(1 + fT) + fT] 
'W ~ 3(1 +fT)^ 


Kp = 3 H^, 


(32a) 

(32b) 


so that (|23 can be related to and expressed as 


y2 = 

which looks like a FLRW Friedmann equation for a source 
made up with dust plus a time dependent A term, p, but 
both given in terms of their averaged densities In 

this context, we can identify the relative velocity Y/Y as 
a sort of averaged Hubble factor. 

















VI. REGULARITY CONDITIONS AND 
ASYMPTOTICS 
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For a spherically symmetric spacetime CHJ character¬ 
ized by the functions R, T, f and F must com¬ 

ply with 

R > 0, (34a) 

1 + fT > 0, (34b) 

l + FT > 0. (34c) 

Hence, from (HU, the condition p >0 becomes 

[i{3TH + f) f F + 6H {f + F)]f > 0 (35) 

so that we can identify t such that R = 0 with a simulta¬ 
neous big-bang singularity associated with p, —>■ oo (with 
p finite) and two singularities for which p ^ oo (with p 
finite) 

iFfT = 0, (36a) 

\+FT = 0, (36b) 

which, in general, will be marked as non-simultaneous 
surfaces in the (t, r) coordinate plane. Notice that 136all 
is a non-simultaneous big-bang, or a central singularity 
(=^ Y = 0), while (I36h|l is a shell crossing singularity 
(=> Y' = 0), both analogous to singularities in Lemaitre- 
Tolman-Bondi and Szekeres dust solutions Hence, 
we will demand that 


Y' > 0 holds for all (t,r) complying with Y > 0, (37) 

or equivalently, that doe s not occur in the coordi¬ 

nate range that complies with 134bll . As shown in figure 
1 , condition 123 can be be satisfied by suitable choices 
of free parameters. 

A usual regularity demand in spherically symmetric 
spacetimes is the existence of a symmetry center, char¬ 
acterized by the regular vanishing of the surface area of 
the orbits of SO (3). This condition defines the worldline 
of a privileged “central” observer for whom the universe 
appears isotropic. From m and marking this worldline 
by r = 0, we must have for all t 

Y(,) = Y(,) = 0 (38) 

where Y(c) = Y(t, 0). We shall use henceforth the 
subindex (c) to denote evaluation at r = 0. Since spacial 
gradients of p, 0, E must vanish at the center, equations 
(nni and imply that /(,,) must be bounded and 

= 0, so that F(c) = /(c). The central values of these 
quantities are then given by 


K p(c) 


4 3g[l + /(c)r] + /(c)r 

[l + /(c)T]2 


(39a) 


Q(c) 

3 

S(c) 


’^1 rr, 2 / 0 )^ 

.^J(c) 3[1 + /(c)T]’ 

0 [cr“ b](c) = 0 


(39b) 

(39c) 



FIG. 1: The plot depicts the coordinate representation of 
the non-simultaneous “big bang” singularity Y (r, r) = 0 as 
a “surface” in the (r, r) plane, with the dimensionless time 
T defined by (55a) and other free parameters given by (54), 
(55) and (60). The domain of regularity, Y > 0, is shown 
with worldlines of fundamental observers (arrows). Notice 
that the surface Y' — 0 occurs in the region with negative Y 
(gray shading). 


where /(c) > 0 and /(r) for r > 0 must be selected so 
that p decreases with increasing r. 

In order to get an idea of the behavior of the models in 
the asymptotic limit t —> oo, we may assume as t —f oo 
an asymptotic power law scaling R ^ for fc > 0. From 
(HHi, we have T ^ 0 in this limit, but T tends to a finite 
asymptotic value (T ^ T*) only for k > 1/3. From 114II . 
113 we have 




3 


k{2- 3k) 


(40) 


while, assuming / and F everywhere finite, we have from 
da and the ratio 


0, for k > 1/3, (41a) 

const, for 0 < fc < 1/3. (41b) 

Therefore, the homogeneous DE fluid dominates asymp¬ 
totically over the DM component even if p > 0 (or equiv¬ 
alently k < 2/3), while values of k associated with nega¬ 
tive p always yield the asymptotic ratio (ETall . It is not 
difficult to verify that the same asymptotic behavior oc¬ 
curs if R follows exponential or logarithmic asymptotic 
scalings, R exp(fc t) or R ^ In t. Thus, it is generi- 
cally possible to have an asymptotically DE dominated 
asymptotic scenario for all fundamental observers, which 
is an important property of the models under consider¬ 
ation. As long as R scales asymptotically faster than 
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R ^2/3 are finite everywhere, we have p < 0 

together with 

^ > iJ, E ^ 0, as t ^ oo, (42) 

so that, regardless of the choice of / (spacial dependence), 
the models homogenize and isotropize for all fundamen¬ 
tal observers {R constant) as t oo, implying that all 
fundamental observers detect in this asymptotic limit lo¬ 
cal conditions that are very close to those of a FLRW 
cosmology characterized by a DE source with fi and p. 

It is also interesting to examine the behavior of the 
models as r —> oo along hypersurfaces of constant t, the 
rest frames of fundamental comoving observers. It is im¬ 
portant to notice that the radial coordinate, r, merely 
labels the fundamental comoving observers, so it has no 
invariant meaning and can be though of as a dimension¬ 
less ratio, such as r = £/£□, where £o is the present value 
of an arbitrary length scale (for example, the scale of ho¬ 
mogeneity £o ~ 100 — 300 Mpc). Thus, we can consider 
“local” scales {£ £o) as marked by r I, so that 

asymptotic “cosmic” scales correspond to r I, while 
the “transition” from local to cosmic scales is roughly 
given by r ^ I. 

Since / is the free function that governs the depen¬ 
dence on r, the variation of all quantities in the radial 
direction is tied to specihc choices of this free function. 
In particular / should be selected so that p decays as r 
increases. Also, the spacial dependence of AI, the “ef¬ 
fective mass” of DM defined by is governed by the 
term r^f, thus it is also convenient to demand that / > 0 
and that r^f must be a monotonously increasing func¬ 
tion. Still, considering all these restrictions, it is possible 
to select / so that 


/ —> 0, F ^ 0, as r ^ oo, (43) 

which is a sufficient condition for 

^ ©|r=cx3 ^ H, p|r=oo ^ 0, E|r=oo ^ 0. (44) 

Thus, since p, is independent of r, these limits imply that 
local conditions for fundamental observers with large r 
are similar those of fundamental observers of a FLRW 
cosmology whose source is a fluid with p and p (that is, 
a “pure” DE fluid). Another possibility is furnished by 
the choice: 


f ^ f*, F ^ /*, as r ^ oo, (45) 
where f* < /(c) is a positive constant. This yields 


^p|r—ooC — ^/ 


2 _ 4 3iL[i + /*r] + /*r 


[l + f*TY 


f, (46a) 


©|r—CX3 — 


= H 


2f*T 


r, (46b) 


= 0 


3 [1 -f /* T] ’ 

[cr“ h]r=oo = 0. (46c) 


The preference of one of these choices of asymptotic be¬ 
havior along the rest frames depends on the problem one 
is interested to study: if we want to examine a large 
scale (supercluster scale or larger) spherical inhomogene¬ 
ity whose evolution requires that we somehow “plug in” 
the effects of a cosmological background, then the choice 
(im may be preferable, while the choice m may be 
preferable for a relatively small scale and/or large den¬ 
sity contrast description of an homogeneity (cluster of 
galaxies) that ignores cosmic effects (see figure 3). 

It is also possible to perform a smooth matching, for 
a comoving radius rj,, of a region 0 < r < of the 
DM and DE mixture to a spatially flat FLRW spacetime 
characterized by p and p, occupying r > Vh- Necessary 
and sufficient conditions follow if /(;,) = F'(j) = 0, where 
(b) denotes evaluation at r = r;,, so that p{p) = E(;,) = 0 
and 0(b)/3 = FI hold for the FLRW region. However, 
such a matching also requires the mass function M in 
to be continuous (at least C^) at r = rb, which from 
implies M(h) = 0, but Ad is a volume integral of p. 
Therefore, if this integral must vanish at r = Cb for all t, 
the dust density p in the integrand must necessarily be 
negative in a Hnite domain of the mixture region 0 < r < 
Tb. Hence, we will not consider this type of matching any 
further. 


VII. OBSERVATIONAL PARAMETERS 


The quantity FI in is the Hubble expansion factor 
associated with a FLRW geometry, for the inhomoge¬ 
neous metric CHI the prop er generalization of this pa¬ 
rameter is given by [ 13 , [ill 

U = ®-H(Tabn“n^ (47) 

where the vector n“ complies with = I, UaTF = 0. 
For a spherically symmetric spacetime, it is necessary to 
evaluate n“ for general comoving observers located in an 
“off-center” position in the spherical coordinates (r, 0 , (/) 
centered at r = 0. For the metric equation EH) 
becomes in general 

n = H+FT 

^ ^ 2 [/(l + Er) + f(f-/) cos^y,] 

3(i + /r)(i + Er) ’ 


where '0 is “observation angle” between the direction of 
a light ray and the “radial” direction for a fundamental 
observer located in (r, 6, ip) El. Therefore, the exact 
local values of the observational parameters D for DE 
and DM are 


n 


DE 


n 


DM 


K /i 

K pc^ 

3T-F 




H + FT 


K pc 


2 > 


H + FT 


pc 

P 


Dr 


(49) 

(50) 
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while the acceleration parameter is na 


q = 




+ 


^DE “t” ^E 


1 + 


3p/m 

1 + p(? I [i. 


(51) 


where S is given by ll^ . 

If we consider the flow of cosmic DM with density p 
at the length scale of the observable universe 3000 h 
Mpc), then present day values of shear and DM density 
gradients in comparable scales are severely restricted by 
the near isotropy of the CMB ^^1 


\\<Jab 


’V 6 |S|' 

0 

0 

0 


< 10 “^ 


(52) 


'ha P,b 


\-] 

. p . 

0 

. p . 


< 10“^ (53) 


Hence the large scale spacial dependence of the observa¬ 
tional parameters lITTl . CT and (ED must also be 
restricted by these bounds. However, these restrictions 
can be strongly relaxed, at a local level, if we examine 
the spacial variation of local values of DM density and 
observational parameters in scales smaller than the ho¬ 
mogeneity scale ~ 100 — 300 Mpc. As mentioned in the 
previous section, scale considerations influence the choice 
of boundary conditions or (B51) . 


VIII. A SIMPLE EXAMPLE: THE “GAMMA 
LAW”. 



In order to illustrate how to work out the expressions 
we have derived and how to calculate relevant quantities 
of the models, we consider now the simple case of a ho¬ 
mogeneous DE fluid satisfying a simple equation of state 
known as the “gamma-law” 

P = ( 7 - (54) 

where 7 is a constant. The dust plus homogeneous 
fluid mixtures that we are studying were examined previ¬ 
ously assuming (among other choices) this equation 
of state, but placing especial emphasis in a dust and ra¬ 
diation (7 = 4/3) mixture. Since our emphasis is now on 
modelling DE sources, we will assume 0 < 7 < 2 / 3 , so 
that —1 < p//i < —1/3. In this case we have from (II 311 . 
(trai . (HU), ((13 and ED 

a = ^ ^ ^ 7/1 Hot, (55a) 

ito ^ 

o ej-2 ^ (^-^ 0 )^ 

K. = 3H = - - -, (55b) 

where Hq — 100km/(secMpc), = 0.7 and T* is a di¬ 
mensionless constant denoting the asymptotic value of T 
(we have then in (1 1511 the choice cq = § 7 / 1 ). 


FIG. 2: Figure (a) illustrates how the present DM density 
profile, p/p,, exhibits an increasingly large inhomogeneous 
contrast as 5 ranges from 5 = 670 (maximal contrast) to 
(5 = 0.1 (minimal contrast). We have used (54) and (55) 
and the free parameters /, 7 , /*, T* are those given by (60) 
and (64). Notice how for all 5, we have p > pc^ in the cos¬ 
mic background region, with all curves tending as r —> 00 to 
the cosmic ratio Ddm/Dde ~ 0.56. Figure (b) displays the 
plot of a measure of density contrast between the center and 
r = 00 , given by Ap/p = p{to, 0)/p{to, 00 ) — 1, as a function 
of logjQ 5. Notice how the contrast grows proportionally to S, 
well beyond the linear perturbation regime. 


Notice that the assumptions (HD have been obtained 
from the FLRW equations m and El and yield a 
power law form for the function a (equivalent to the 
FLRW scale factor). Therefore, following this form 
of the homogeneous DE fluid is equivalent to a scalar 
field with an exponential potential following the so-called 
“scaling law”. 

The density of the dust component and the generalized 
Hubble factor are found by inserting (HD into (13 and 

(33 
































FIG. 3: The panels depict the logarithmic profile of p/p(^c), 
with = p(r, 0 ) for a sequence of values of constant r, up 
to the present t = tq. The free parameters 7 , 5, T* are those 
given by (64). In panel (a) we selected /* = 100, hence p/Pc 
tends to the nonzero asymptotic cosmic DM density given by 
(46a). We see a power law decay of DM density only in the 
transition (r ~ 1 ) from the central overdensity region to the 
cosmic background. In panel (b) we selected /* = 0, thus the 
curves closely resemble isothermal profiles, in which p decays 
to zero approximately as r“^. This is not surprising since 
A4 —> r asymtotically with the choice (60). 


Kpc^ 


U 


3 7^ Hi [2 T* / f + (1/7) (/ + F)] +j,fF 

t 2 [(i + /T*)ri/7i _^i/] [(i + i^r*)Ti/7i 

hHo _ [(l + FT*)ri/7i / + |(F-/)ri/7i cos^V'l 

r I ^ [(i + /T*)t1/7i _^i/] [(l + i^r*)ri/7i j ■ 


2-37 

2-7 ’ 


(56) 

(57) 


From ( 1 ^ . we see that R scales as so that fc > 1/3 
for 7 < 2 , hence for the 7 values that we are interested 
we should obtain the ratio pc? jp given by (I41all . Using 
(I55bll and (1^ and assuming an arbitrary but finite / 


and F, we obtain in the limit r ;g> 1 


p? 2 r*/f+ ( 1 / 7 ) (/ + u) 37 ^ 
(i + /r*)(i + FT*) ri /71 
hHo 


0, (58) 


n 


T 


= H 


(59) 
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(a) (b) (c) 



FIG. 4: The asymptotic value of / in (60) is given by /*. Assuming T* — 0 and an arbitrary S, we can find a suitable value 
for /*, for any given 7 in (54), by demanding that the observational parameters Ude, S4dm, q (panels (a),(b) and (c)) have 
appropriate “present” cosmological values in the ranges (63) as r —> 00 for t = tq (gray stripes). Each curve is marked by a 
given value of 7 . Notice that for 7 = 0.15, the choice f* ~ 100 yields Ude = 0.64, Udm = 0.35, q = —1/2, which are reasonably 
close to currently accepted observational data. For 7 closer to 0 (cosmological constant), we would have to select larger values 
for /*, while larger 7 close to 0.3 correspond to f* ~ 1. 


indicating that for all cosmic observers the mixture ho¬ 
mogenizes and isotropizes as the homogeneous DE fluid 
dominates asymptotically over the cold DM component. 


IX. NUMERICAL EXPLORATION. 


Having found T, H and a = R/Rq for the particular 
case of a gamma law we only need to select the 
function / = /(r) in order to render the models fully 
determinate. A convenient form for / is 

/ = r + T^. ( 60 ) 

1 -I- r"' 

so that both type of asymptotic boundary conditions, 
or 63-1111), can be accomodated by selecting 
/* to be zero or nonzero. Notice also that asymptotically 
as r ^ 00 , we have: A4 ^ r if f* = 0 and A4 ^ if 

r >0. 

Considering 63, the regular evolution range for the 
models is the coordinate range (r, r) where 63 holds 
with a, H, T given by 63,1^ and /, F given by 63 
and unj). Thus, assuming as the initial “big-bang” the 
coordinate surface (I36all marking Y = 0, we can take 
the big bang time as the value r = run at this surface 
corresponding to r ^ 00 (see figure 1 ), that is: 


7i (/* + </) 

l + {f*+S)T* 


(61) 


Thus, considering the “age of the universe” roughly as 
Ato « 14 Gys and h « 0.7, the present cosmic era corre¬ 
sponds to 


To = Tbb + - jhHo Ato « Tbb + 3.17 7 


(62) 


As shown in figure 1 , the big bang surface 136all is not 
simultaneous, thus for any hypersurface r = constant, 
the regions near the center at r = 0 will be “younger” 
than those asymptotically far at large values of r. 

The parameter S = f(^c) — f* > 0 in (11)01) provides the 
a measure of inhomogeneity contrast, or “spacial” varia¬ 
tion of all quantities along the rest frames (t = constant) 
between the symmetry center r = 0 and r 00 . Thus, 
a sufficiently large/small value of <5 makes the values at 
r = 0 and r 00 sufficiently close/far to each other, 
thus indicating small/large “contrast” or degree of inho¬ 
mogeneity. We test the effect of 6 along the present day 
hypersurface r = tq by plotting in figure 2 a the spacial 
profiles of pc^//r for a sequence of <5 values, while figure 
2 b depicts the density contrast between the overdensity 
at center and the cosmological background, as a function 
of logio 

The effect of selecting /* zero or nonzero is shown in 
figures 3a and 3b, by means of a comparison between 
the resulting DM density profiles that follow by evalu¬ 
ating 63 and 63 along a sequence of hypersurfaces of 
constant r (up to tq). The asymptotic condition ( 11 !^ 
(/* = 0) yields an approximately r~^ power law decay 
of p that resembles a standard isothermal profile, all of 
which is consistent with (63 and with the fact that the 
form of / in 63 , /* = 0 , leads to A4 —> r, the 

same asymptotic behavior of the standard “isothermal 
sphere”. On the other hand, the choice ll^ (/* > 0 ) 
leads also to a power law decay of p, but towards a cos¬ 
mological background with asymptotic density given by 

63 - 

The appropriate numerical value for the asymptotic 
constant, /*, can be found by demanding that the cos¬ 
mological observational parameters Ode, 0 dm, 9 , evalu¬ 
ated in the cosmic background (r —> 00 ) at the present 
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FIG. 5: The figure illustrates the effects of anisotropy 

through the dependence of the profiles of Ude and Udm, along 
T = To, on the off-center “observation” angle ')/>. Free param¬ 
eters are given by (64). The thick black curves correspond to 
■(/> = 0 , 7 r, while the gray curves provide the profiles of Ude 
and SIdm for various other values of ip- In particular, the 
larger “peaks” of these parameters occur for ip = 7 r/ 2 , 37 r /2 
and around r ~ 1. Notice how these effects of anisotropy are 
negligible in the overdensity region around r = 0 and in the 
cosmic background for large r. For the deceleration parameter 
q these effects are negligible for all r. 

era (r = tq), take (for a given 7 ) reasonably close val¬ 
ues to those currently accepted from observational data. 
Since (1^ with T* = 0 and / given by (IHHll . evaluated 
at T = To and r ^ 00 , is independent of 8 and ip, we can 
plot floE, riDMj <1 as functions of f* and 7 . As figures 4a, 
4b and 4c illustrate, the desired value of /* for any given 
7 can be selected so that the forms (03, (EDI and lEU 
at r = To and r —> 00 yield: 

0.6 < Ode < 0.7, (63a) 

0.3 < floM < 0.4, (63b) 

-0.5 < q< -0.4 (63c) 

In particular, if we select 7 = 0.15, an appropriate value 
is f* « 100, leading to Ode ~ 0.64, 0dm ~ 0.35 and 
q « —0.5. Notice that Ode + 0dm ~ 1, but the present 
Omega for DM would be slightly higher than the cur¬ 
rently accepted value 0 dm ~ 0.3. Thus, we could argue 
that these parameter values would become a very accu¬ 
rate approximation to actually inferred cosmological pa¬ 
rameters if we would consider p as the compound density 
of DM and baryonic matter. 

The efects of anisotropy emerge in the dependence of 
TL, as given by 03 , on the off-center “observation an¬ 
gle” Ip. This implies dependence on ip for Ode, 0dm, 9 - 
Considering the free parameter values 

T* = 0, f* = l00, 5 = 200, 7 = 0.15, (64) 

complying with the cosmic background ranges 103 , fig¬ 
ure 5 displays Ode and 0 dm, evaluated at r = rg, as a 


FIG. 6: The figure illustrates the effects of inhomogeneity by 
showing how the time evolution of Ode and 0 dm (for ip = 0) 
is very different along the center of the overdensity at r = 0 
(crosses) and in the cosmic background (r —» 00 , solid curves). 
The free parameters are given by (64). Notice how in the 
present cosmic era (r = tq) DE dominates over DM in the 
cosmic background, but DM dominates over DE in the center. 


function of log^g r for assorted fixed values of ip. The 
same profiles of Ode and 0dm occur for ip = 0 and 
Ip = TT (thick black lines) and, in general for any two 
values of ip that differ by a phase of tt, with the high¬ 
est “peaks” corresponding to gray curves with ip = 7r/2 
and Ip = (3/2)7r. This singles out two “preferential” dis¬ 
tinctive directions: one along the axis ip = 0, tt and the 
other along ip = 7r/2, (3/2)7r. This is a clear representa¬ 
tion of a quadrupole pattern, as expected for a geodesic 
but shearing 4-velocity |l3,[i3- Also, as revealed by fig¬ 
ure 5, the curves for the various ip differ from each other 
only in the transition region near r ^ 1, thus the effects 
of this quadrupole anisotropy are negligible near the cen¬ 
ter of the local overdensity and in the cosmic background 
asymptotic region. 

The effects of inhomogeneity are illustrated by figure 6, 
for the same parameters lIHUl but keeeping ip = 0 fixed, 
by plotting Ode and Udm as a function of r for r = 0 
(crosses) and r = 00 (solid lines). Notice how for the 
present cosmic era, r = tq, we have Ode dominating at 
the cosmic background region (r —> 00 ), but IIdm dom¬ 
inates in the overdensity region (r = 0). This effect of 
inhomogeneity can be further appretiated in figures 7a, 
7b and 7c, displaying Ode, 0dm and q as functions of r 
and logj^Q r for ip = 0, tt and using the free parameters 
03 . It is particularly interesting to remark how as r 
grows we have: Ode ^ 1,0dm ^ 0 and q —1/2 
for all r, as expected for a DE dominated asymptotically 
future scenario and associated with an ever accelerating 
universe that follows a “repulsive” dynamics. However, 
in the present cosmic era (thick black curve) this repul¬ 
sive accelerated dynamics on which Ode dominates and 
9 < 0 only happens in the cosmic background region, 
with 0dm > Ode and 9 > 0 { i . e . “attractive” dynamics) 
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FIG. 7: The figures illustrate the full dependence of 

SIde, Udm and q on t and r (panels (a), (b) and (c)) for 
Ip = 0 and with free parameters given by (64). The thick 
black curves denote the hypersurface of present cosmic time 
T = To. The figures clearly show the differences among the 
overdense region around r = 0, the cosmic background for 
large r and the transition zone between them around r ~ 1. 
As in the previous figure, it is evident that DE dominates over 
DM in the cosmic background, but DM dominates over DE 
in the center, with both converging as r ^ oo to asymptot¬ 
ically homogeneous states with values Ude ^ 1, Ddm ^ 0. 
Notice in (c) how the deceleration parameter, q, is negative in 
the cosmic background (accelerating universe), but is positive 
in the overdense region where the dynamics of local gravity 
should not be repulsive. 


in the local overdensity region with a relatively large DM 
density contrast 6 = 200. 

Conditions (EH) and E3 place stringent limits on large 
scale deviations from homogeneity and anisotropy, but 
these bounds do not apply to local values of these quanti¬ 
ties. It is interesting to plot these quantities for the cases 



FIG. 8: Level curves of the ratio \/6S/0. Free parameters 
are those given by (64). Notice that condition (52) is satisfied 
for most of the range of r at r = tq, failing to hold only in 
the transition (r ~ 1) between the overdensity and the cosmic 
background. 



FIG. 9: The interaction function Q(r,r) given by (66). Free 
parameters are given by (64) and the level curve Q = 0 is 
the thick black curve. Notice how Q > 0 in the present cos¬ 
mic era (DE transfers energy to DM), then changes sign (DM 
transfers energy to DE), with Q —> 0 for large r, thus indicat¬ 
ing an asymptotic state in which DE and DM only interact 
gravitationally. 


depicted in the previous figures, all characterized by ED), 
(iCTi . (1^ and (IHUl . Figure 8 depicts level curves of the 
logarithm of the shear to expansion ratio of E3), with 
0 and S given by (ETll-(l^. Even considering the rela¬ 
tively large value 6 = 200, condition ESi holds througout 
most of the coordinate range (r, r) including the far range 
“cosmological background” region of large r and the lo¬ 
cal overdensity region near the symmetry center r = 0, 
so that for the present cosmic time r = tq it only ex¬ 
cludes the relatively small scale local region around r ~ 1 
that marks the “transition” from the local overdensity to 
the cosmic background. However, a choice like 6 = 0.01 
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FIG. 10: Qualitative logarithmic plot of the coupling function 
defined by (71), with ^ oc 1/r^ and free parameters 
given by (64). The hypersurface r = tq is displayed as a thick 
black curve. Notice the clear power law time dependence 
of just slightly modulated by the spacial dependence. 

This indicates that x = x(4>) holds approximately for a scalar 
field associated with (54) and (67). 


would yield similar level curves, but with values three 
orders of magnitude smaller, thus denoting a state of al¬ 
most global homogeneity, since (PI would hold in almost 
all local scales in r = tq. A graph that is qualitatively 
very similar to that of figure 8 emerges for condition (15811 . 

This difference between the dynamics of local inho¬ 
mogeneities and that of the cosmic background cannot 
be appretiated in such a striking and spectacular way if 
one examines DE and DM sources by means of the usual 
FLRW models and their linear perturbations. 


X. INTERACTION BETWEEN THE MIXTURE 
COMPONENTS 


As we mentioned before, the two mixture components: 
DM (inhomogeneous dust) plus DE (homogeneous fluid) 
are not separately conserved. Considering m and (EH, 
the energy balance for the total energy-momentum ten¬ 
sor, e -I- (e -I- p) 0 = 0, can be written in terms of the 
DM and DE components as 


p-hp(3i?-kZr)] -b [/i-b (/i-kp) (siT-bzr) 


= 0 , 


. 

Since each term in square brackets in the left hand side of 
(EHl corresponds to the energy balance of each mixture 
component alone, a self-consistent form for describing 
the interaction between the latter is given by 


fi+in + p) 

P + P 


3H + ZT 
+ Zf 


= -Q, 
= Q, 


(66a) 

(66b) 


where Q = Q{t,r) is the interaction term. Since the 
physics behind DM and DE remains so far unknown, we 


cannot rule out the existence of such interaction. No¬ 
tice that once a given model has been determined by 
specifying an “equation of state” p = p(p,) and a form 
for / = /(r), as we have done in the previous sections, 
this interaction term would also be fully determined. In 
general, if the interaction term in (TO is a negative val¬ 
ued function, then DM transfers energy into the DE and 
viceverza. Considering the free parameters given by iE3J), 
(PI and (IMI) . we plot in figure 9 the interaction term Q 
in (ra . as a function of r and log^Q r. Notice that Q is 
initially positive and remains so today (DE trasfers en¬ 
ergy to DM at r = tq) but will change sign in a future 
time (DM trasfers energy to DE), tending to zero asymp¬ 
totically as r ^ oo. Thus, the time-asymptotic state is 
that of only gravitational interaction between DE and 
DM {i.e. separate conservation of each component). 

However, the relevant question is not so much the ex¬ 
plicit functional form of Q, but its interpretation in terms 
of a self-consistent physical theory that would be regulat¬ 
ing the interaction between DM and DE. In fact, one of 
the challenges of modern cosmology is to propose such 
a self-consistent theoretical model of this interaction, 
while agreeing at the same time with the experimental 
and observational data. In this context, the interaction 
between DM and DE has been considered, using homo¬ 
geneous FLRW cosmologies, in trying to understand the 
so-called “coincidence problem”, that is, the suspiciously 
coincidental fact the DE and DM energy densities are of 
the same order of magnitude in our present cosmic era 

[I i 11 ii 111 113 

If the homogeneous DE fluid corresponds to a 
quintessense scalar field, (p = Pit), with self-interaction 
potential V{(j)), we have instead of II54II : 


P=^ + V{P), p=^-V{4>). (67) 

In this case, the interaction term in (TO) can be associated 
with a well motivated non-minimal coupling to the DM 
component. Consider a scalar-tensor theory of gravity, 
where the matter degrees of freedom and the scalar field 
are coupled in the action through the scalar-tensor metric 


^ST = I rf"a:v1^{|-i(V<^)2 

+ ( 68 ) 

where xi4>)~'^ is tbe coupling function. Cm is the matter 
Lagrangian and v is the collective name for the matter 
degrees of freedom. Equations become 


P + P 


SH + ZT 
3H + Zf 


_P_X 

# 20 x’ 

P X 

'2 x' 


(69a) 

(69b) 


so that, the coupling function x('/') the interaction 
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term Q are related by 


Q = -f| (70) 

Therefore, once we determine a given model, so that Q 
can be explicitly computed, we can use m to find the 
the coupling function y that allows us to relate the un¬ 
derlying interaction with the theoretical framework asso¬ 
ciated with the action 

For the models under consideration, we can integrate 
(in general) the constraint 17011 with the help of dia, 
(P|l. IIH), (EIJ, 113 and using 0 = {d/dt) [ln{Y‘^ Y')]. 
This yields 

y-i/2 = ar)pY^Y' 

(71) 

where ^(r) is an arbitrary function that emerges as a 
constant of integration. Notice that the models require 
4> = 4>(t), so that the assumption y = xi'P) implies 
X = x(0- However, from 17111 . we have (in general) 
y = y(t,r), with the case y = y(t) occurring for the 
following particular cases, associated with very special 
forms of / and ^ 


F = 

0 , 

/ocr-3/2, 



x“ 

oc H, 

(72a) 

f + F = 

0 , 

f oc r“3, ^ oc 



x“ 

oc f + 3TH 

(72b) 


or, if F 7 ^ 0 and f+F ^ 0, then / and ^ must be obtained 
from the constraints: 


3/" /'" 

2 / + /2 


5 = 0 , 

r 


i: 


f 

6/ + 6rr+ r" / 

r{r f' + 3/) 


(73) 


However, since m and 123 yield very special forms of 
/, ^ and of y, we prefer to apply under the most 
general assumption that the coupling function y should 
be a function of (j) and of position, i.e. y = y( 0 (t),r), 
as given by GD for suitable forms of the free functions 
/, H and T (thus, (j) and V{(j))), hence ^ can be consid¬ 
ered a wholly arbitrary function. Considering the free 
parameters given by 113, itSBl) . (Uni) and (IMIl . and tak¬ 
ing ^ oc 1 /r^, figure 10 illustrates qualitatively how the 
coupling function y“^/^ given by G3 follows a power 
law time dependence, slightly modulated by the change 
of local DM densities from the local overdensity region 
to the low density cosmic background. Since the scalar 
field (IHTIl that corresponds to the equation of state ll?)^ 
and the parameters also exhibits a power law time 


dependence M , we do have approximately y = as 
required by ra . 

However, figure 10 is just a qualitative plot. We should 
point out that relating the interaction term, Q, to the for¬ 
malism represented by is strictly based on the formal 
similitude between the field equations derived from the 
action P|). on the one hand, and equations and 
on the other. Also, the interaction between DM 
and DE in the context of is severe ly c onstrained by 
experimental tests in the solar system A more de¬ 
tailed and carefull examination of the relation between y 
and (j) that incorporates properly these points, as well as 
the application of to the models presented here, will 
be undertaken in future papers. 


XI. CONCLUSION. 

We have presented a class of inhomogeneous cosmo¬ 
logical models whose source is an interacting mixture of 
DM (dust) and a generic DE fluid. The relevance of 
the present paper emerges when we realize that there 
are surprisingly few studies in which DE and DM are 
the sources of inhomogeneous and anisotropic spacetimes 
(see ll^), as practically all study of the dynamical evo¬ 
lution of these components is carried on in the context 
of homogeneous and isotropic FLRW cosmologies or lin¬ 
ear perturbations on a FLRW background. There are 
also very few papers that examine the possibility of non- 
gravitational interaction between DE and DM. 

Once we assume or prescribe an “equation of state” 
(i.e. a relation between pressure, p, and matter-energy 
density, p, of the DE fluid), we have a specific DE model 
(quintessense scalar fields, tachyonic fluid, etc) and all 
the time-dependent parameters can be determined by 
solving two coupled differential equations reminissent of 
FLRW fluids. Since the spacial dependence of all quanti¬ 
ties is governed by the function / = /(r), once the latter 
is selected the models become fully determined (though, 
various important regularity conditions must be also sat¬ 
isfied: see section VI). In order to work out this pro¬ 
cess we chose the simple “gamma law” equation of state: 
p = (7 — l)fi (equation iHtHi L leading to analytic forms 
for all relevant quantities, including the main observa¬ 
tional parameters, Hde, and q. Our choice for a DE 
fluid complying with is equivalent to a scalar field 
with exponential pote ntial, satisfying a scaling power law 
dependence on t [T^. Although this is a very idealized 
quintessense model, our aim has been to use it as a guide¬ 
line that illustrates how more sophisticated DE scenarios 
can be incorporated in future work involving the models. 

As we have mentioned in section VI, the models homo- 
geneize and isotropize asymptotically in cosmic time for 
all fundamental observers and/or assumptions on the DE 
fluid, thus they are well suited for studying the interac¬ 
tion between DE and DM in the context of the evolution 
of large scale inhomogeneities (of the order of the scale of 
homogeneity ~ 100 —SOOMpc). By selecting appropriate 
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boundary conditions, we can examine inhomogeneities at 
various scales and/or asymptotic conditions (see figures 
2 and 3). In particular, we have explored the case of a lo¬ 
cal DM overdense region, whose scale can be arbitrarily 
fixed and with an asymptotic behavior that accurately 
converges to a cosmic background characterized by ob¬ 
servational parameters that fit currently accepted obser¬ 
vational constraints: 0.6 < Dde < 0.8, 0.2 < Ddm < 0.4 
and q « —0.5 (see figure 4). As illustrated by the various 
graphical examples that we have presented, this inter¬ 
play between a local overdensity, a cosmic background 
and a transition region between them, shows in a spec¬ 
tacular manner how inhomogeneity and anisotropy lead 
to interesting and important information that cannot be 
appretiated in models based on FLRW metrics and/or 
linear perturbations. For example, as revealed by figure 
5, the effect of anisotropy emerges as a dependence of 
observational parameters on local observation angles in 
“off center postions”, an effect which is only significant 
in the transition between the overdensity and the cos¬ 
mic background. On the other hand, as shown by figures 
6-8, inhomogeneity allows for radically different ratios 
between DM and DE in the overdensity and the cosmic 
background, so that DM dominates over DE, locally, in 
the overdense region, as a contrast with DE dominat¬ 
ing DM, asymptotically, in the cosmic background (as 
expected). Also, while q is negative in the cosmic (DE 
dominated) background, thus denoting the expected “re¬ 
pulsive” accelerated expansion at large scales, we have 
q > 0 along smaller scales in the local overdensity. For 
all parameters there is a smooth convergence between 
local and asymptotic values in the transition region. 

We have also examined the non-gravitational interac¬ 
tion between DE and DM. Ploting the interaction term, 
Q, (figures 9) shows that energy flows from DE to DM 
at the present cosmic era, with the flow reversing direc¬ 
tion in the future and evolving towards an asymptotic 
future state characterized by pure gravitational interac¬ 
tion: Q ^ 0. If we take the DE fluid to be a quintessense 
scalar field, the DE vs DM interaction can be incorpo¬ 


rated to the theoretical framework of an action like ra . 
associated with a non-minimal coupling of scalar fields 
and DM. Since DM is inhomogeneous while the scalar 
field is homogeneous, only for some particular forms of 
spacial dependence (ie. the function /) we obtain a 
coupling function expressible as y = xi4>)- In general, 
we have to allow for the possibility that y = x((/)(t),r), 
but as shown by figure 10, the spacial dependence of y 
under the assumption 15411 (scalar field with exponential 
potential) mantains a power law dependence that is qual¬ 
itatively very similar to that shown in the homogeneous 
case. However, we have examined this interaction just 
in qualitative terms, with the purpose of illustrating the 
methodology to follow in future applications. 

As guidelines for future work, we have the application 
of the models to more sophisticated and better motivated 
DE formalisms, perh^s in the context of the “coinci¬ 
dence” problem [a, IE 0 IS @> Q • An interesting de¬ 
velopment would be the study of the case AT 0, in which 
hypersurfaces of constant t (rest frames) have nonzero 
curvature (see section II). This requires solving the non¬ 
linear equation ©, which would certainly need numerical 
techniques and perhaps assuming an homothetic symme¬ 
try that could reduce it to an ordinary differential equa¬ 
tion. The advantage would be the possibility of describ¬ 
ing the effects of DE on structure formation, since a DM 
overdense region with K > 0 could collapse, locally, as a 
DE dominated cosmic background with K < 0 expands. 
Finally, we have examined the DE/DM mixtures that 
can be constructed with the Szekeres-Szafron models of 
“class 1”, as those of ^3, thus it still remains pending 
the study of “class 2” models like those of |T^. 
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